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A recent investigation of the single spin asymmetry (SSA) in low virtuality electropro- 
duction/photoproduction of J/ip in color evaporation model is presented. It is shown 
that the asymmetry is sizable and can be used as a probe for the still unknown gluon 
Sivers function. 
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1. Introduction 

Large single spin asymmetries (SSA) observed when an unpolarized beam of elec- 
trons or protons is scattered off a transversely polarized target can be explained 
with the inclusion of k± dependence in parton distribution functions (pdf's) and 
fragmentation functions(ff's) ^ l n recent years, there has been a lot of interest in 
investigations of transverse single spin asymmetries in high energy QCD processes 
as they provide information about spin and orbital angular momentum carried by 
the quarks and gluons. Large SSA's have been measured in pion production at 
Fermilab ^ as well as at BNL-RHIC inpp^ collisions^. SSA's have also been ob- 
served by HERMES 4 and COMPASS EI collaborations, in polarized semi-inclusive 
deep inelastic scattering. The magnitude of the observed asymmetries have been 
found to be larger than what is expected from perturbative QCD ® It is possi- 
ble to explain this large asymmetry in terms of transverse momentum depenent 
(TMD) parton distribution func tion s. One is led to a generalized factorization for- 
mula called TMD factorization which in some processes has been proved at 
leading twist and leading order ^ and has been argued to hold at all orders. One 
of the TMD functions is the Sivers function which describes the probability of find- 
ing an unpolarized parton inside a transversely polarized hadron. Here, we propose 
charmonium production as a probe to investigate the Sivers function and estimate 
SSA in photoproduction (low virtuality electroproduction) of charmonium in scat- 
tering of electrons off transversely polarized protons. At leading order (LO), this 
receives contribution only from a single partonic subprocess jg — ¥ cc . Hence, SSA 
in e + p' — > e + J/ip + A, if observed, can be used as a clean probe of gluon Sivers 
function. In addition, charmonium production mechanism can also have implica- 
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tions for this SSA and therefore, its study can help to understand the production 
mechanism for charmonium. 

2. Sivers Asymmetry in J/ift Production : Formulas 

There are several models for charmonium production. In the color singlet model 021 
the cross section for charmonium production is factorized into a short distance part 
for cc pair production calculable in perturbation theory and a non-perturbative ma- 
trix element for the formation of a bound state, which is produced in a color singlet 
state. A more recent model of charmonium production is the color octet model 
This is based on a factorization approach in non-relativistic QCD (NRQCD) and it 
allows cc pairs to be produced in color octet states. Here again, one requires knowl- 
edge of the nonperturbative colour octet matrix elements, which are determined 
through fits to the data on charmonium production. On th e other hand, in color 
evaporation model (CEM), which we use in this 

work ,[12][13] a 

statistical treatment 

of color is made and the probability of finding a specific quarkonium state is as- 
sumed to be independent of the color of heavy quark pair. We have used Weizsacker- 
William s (WW ) equivalent photon approximation for the photon distribution of the 
electron t calculate the cross section for the process e+p^ — > e + J/i> + X at 
low virtuality of the photon. The underlying partonic process at LO is "fg — t cc and 
therefore, the only k± dependent pdf appearing is the gluon Sivers function. For a 
complete calculation of photoproduction of J /ij) one has to consider higher order 
contributions and also the resolved photon contributions ^4 Also the gauge links 
or Wilson lines present in the TMD distributions are important at higher order 

According to CEM, the cross section for charmonium production is proportional 
to the rate of production of cc pair integrated over the mass range 2m c to 2mjj 

a = - \ dM—^r 1 

where m c is the charm quark mass and 2m d is the DD threshold, M 2 is the squared 
invariant mass of the cc pair. The WW distribution function of the photon in the 
electron given by 

f y /e{r,E) = -{ 

7r r 







r 


In 






(in*- 







\r 2 ) 


+ 1 


\ m 




+ 2 









< 2 -"Vi=^| } . (2) 



2r \2-r 

where r is the energy fraction of the electron carried by the photon. 

The single spin asymmetry for semi- inclusive process + B — > C + X is defined 

as 

dcA — da^ . . 

An = dot + dVi (3) 

We assume a generalization of CEM expression by taking into account the trans- 
verse momentum dependence of the Weizsacker- Williams (WW) function and gluon 
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distribution function. The numerator of the SSA can be written as 



dydM 2 d 2 q T dydM 2 d 2 q T 



- J [d?k± 1 d 2 k±g]A N f g/p t (x g , fcx g )/ 7 / e (x 7 , fe_L 7 ) 
x 5 2 {k^ + k ±g - q T )a^ cE (M 2 ) (4) 



where y is the rapidity and qr in the transverse momentum of the charmonium; 
A N fg/pt { x g> fc-Lff) is the gluon Sivers function, / 7 / e (x 7 , fcj_ 7 ) is the photon distri- 
bution of the electron, given in the WW approximation. The denominator would 
have a similar expression involving the unpolarized gluon distribution of the pro- 
ton; fg/p(x g , k± g ), for which we use a gaussian form of A:j_ distribution and a similar 
gaussian form for the transverse momentum dependence of the WW function. To 
extract the asym metry produced by the Sivers function, we calculate the weighted 
asymmetry 



sin(0, T -0 s ) _ f d(j> qT [da t - do 4 ] sm(<p qT - </>g) 
A ~ f d(b qT [dat + da±] 



where 4> qT and <ps are the azimuthal angles of the J/ip and proton spin respec- 
tively. For the gluon Sivers function we have used a model in our analysis, which has 
been used i n th e literature to calculate SSA in semi-inclusive deep inelastic scatter- 
ing (SIDIS) ^ anc j Drell-Yan (DY) process ^ (see ^ for details) . The parameters 

IOQI 

are taken from . Other parameters we use are 

(klg) = (kl 1 )=0.25GeV 2 . 



3. Numerical results 

In this section, we present numerical estimates of the Sivers asymmetry in charmo- 
nium production in different experiments. Model I refers to the parametrization in 

and (a) refers to the parametrization of t he g luon Sivers function in terms of an 
average of the u and d quark Sivers function L ^ 1 . The estimates are obtained using 
GRV98LO for gluon distribution functions. We point out that th e scale evolution of 
the TMD's includin g th e Sivers function has been worked out in l^4 | 25 | anc | recently 
it has been noted in <^ that in SIDIS the evolution indeed affects the Sivers asym- 
metry. However, in the model we consider for charmonium production, namely the 
color evaporation model, the only relevant scale is M 2 which is the invariant mass of 
the heavy quark pair. This is integrated between a narrow region, from Am 2 to 4m|j 
irrespective of the center-of-mass energy of the experiment. So the scale evolution 
of the TMDs is not expected to affect the asymmetry too much. 

In Figs. 1 and 2 we have shown the asymmetry (A s ^ i ^ qT < ^ >s ' ) ) as a function 
of rapidity y and q T respectively for COMPASS (^s = 17.33 GeV) and eRHIC 
(y/s = 31.6 GeV) energies. The plots are for model I (a) and (b) for two choices of 
the gluon Sivers function. We obtain sizable asymmetry in the kinematical regions 
of all the experiments for model I (b). The asymmetry is smaller in model I (a). 
For COMPASS as well as for eRHIC the asymmetry increases with y, reaches a 
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Fig. 1. The single spin asymmetry J 4^ n(< * ,T </>s) for the e + p 1 ^ -> e + J/tp + X at COMPASS as a 
function of y (left panel) and qT (right panel). The plots are for model I with two parameterizations 
(a) (solid line) and (b) (dashed line). The integration ranges are (0 < q-r < 1) GeV and (0 < y < 1). 
The results are given at y/s = 17.33 GeV. 
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Fig. 2. 
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The single spin asymmetry A N ^ qT ' 



for the e + p T — > e + J/ip + X at eRHIC as a 
function of y (left panel) and qj- (right panel) . The plots are for model I with two parameterizations 
(a) (solid line) and (b) (dashed line). The integration ranges are (0 < < 1) GeV and (0 < y < 1). 
The results are given at ^/s = 31.6 GeV. 



maximum and then decreases. This maximum is reached at y ~ 0.6 for COMPASS 
and at y « 1.2 for eRHIC for model I (b). The asymmetry increases with qx for 
both models, and for higher values of qr ~ 0.6 — 0.7 GeV, it becomes relatively 
steady. 

In Fig. 3 we have shown a comparison of the y and qr dependence of the asym- 
metry at different experiments, namely at JLab, HERMES, COMPASS and eRHIC. 
Different experiments cover different kinematical regions, and the plots clearly show 
that the asymmetry is sizable, and that SSA's in charmonium production is a useful 
tool to extract information on the gluon Sivers function. 
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for the e + p^^e + J/^ + Xasa. function of 
y (left panel) and gr (right panel). The plots are for model I with parameterization (a) compared 
for JLab (y/s = 4.7 GeV) (solid line), HERMES (y/s = 7.2 GeV) (dashed line), COMPASS 
= 17.33 GeV) (dotted line), eRHIC-1 (y/s = 31.6 GeV) (long dashed line) and eRHIC-2 
(^/J= 158.1 GeV) (dot-dashed black line). 
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